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ABSTRACT 


A  computer  program  haB  been  developed  for  the  analysis  of  the  com¬ 
ponents  of  a  gamma  pulse -height  distribution  composed  of  a  mixture  of 
two  or  more  gamma -emitting  radic  nuclides.  The  program  consists  of  (l) 
identifying  the  largest  gaama  photopeak  and  matching  its  energy  with 
the  largest  photopeak  of  each  of  the  individual  radionuclides  that 
compose  the  mixture,  (2)  subtracting  enough  of  the  individual  radio¬ 
nuclide  gamma  pulse-height  distribution  from  the  mixture  so  that  the 
photopeak  in  question  is  reduced  by  a  certain  small  percentage,  and 
(3)  continuing  the  incremental  subtraction  process  until  the  photo¬ 
peaks  in  the  mixture  can  no  longer  be  identified.  This  computer  pro¬ 
gram  has  several  advantages  over  other  computer  programs  that  have 
been  developed  to  do  the  equivalent  of  a  radiochemical  analysis  on  a 
mixture  of  radionuclides.  It  will  be  especially  useful  in  analyzing 
the  large  number  of  low  activity  mixtures  of  radionuclides  that  are  to 
be  made  in  the  program  to  determine  the  physical-chemical  species  of 
the  radionuclide  debris  produced  by  underwater  nuclear  explosions. 
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SUMMARY 


Problem 

In  a  study  of  the  physical  and,  chemical  species  of  radioactive 
debris  resulting  from  an  underwater  nuclear  explosion,  complex  mixtures 
of  radionuclides,  sometimes  of  fairly  low  activity,  are  to  be  analyzed. 
A  computer  program  for  resolving  the  gamma  pulse-height  distribution 
of  a  weakly-active  mixture  of  gamma-emitting  radionuclides  into  its 
components  is  required. 

Findings 

A  computer  program  was  developed,  and  then  tested  successfully  on 
the  ganmia  pulse-height  distribution  of  a  synthetic  mixture  of  3  radio¬ 
nuclides,  each  having  several  photopeaks  over  a  2  MEV  range.  The 
program  is  based  on  an  Iterative  subtraction  technique  and  is  promis¬ 
ing  because  errors  in  the  analysis  were  small  (0.5  -  2.0  percent), 
despite  the  fact  that  many  possible  refinements  in  the  program  have 
not  yet  been  made. 
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INTRODUCTION 

This  work  Is  part  of  a  program  to  determine  the  chemical  and 
physical  species  of  radionuclide  debris  formed  from  an  underwater 
nuclear  explosion.  This  underwater  nuclear  explosion  is  simulated  in 
the  laboratory  by  exposure  of  a  small  bead  of  uranium  metal  suspended 
in  seawater  to  a  neutron  pulse  from  a  TRXOA  reactor.  The  absorbed 
neutrons  produce  enough  fissions  (~  2  x  1012  in  10  mg  of  U)  In  a  few 
milliseconds  to  vaporize  the  uranium  and  produce  conditions  that  are 
similar  to  those  of  the  actual  underwater  event. 

Study  of  the  long-lived  radionuclide  debris  requires  making  many 

separations  end  Identifications  on  small  aliquots  of  a  solution  that 
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contains  a  total  of  «.  2  x  10  fissions  In  five  ml  of  seawater.  Radio¬ 
chemical  separation  and  Identification  would  be  time  consuming  and 
difficult  to  do  sequentially  (which  would  be  necessary  because  of  the 
low  activity)  for  several  radionuclides  In  one  sample. 

A  possible  alternative  to  making  these  difficult  radiochemical 
separations  would  be  to  spike  the  uranium  bead  before  Irradiation  with 
relatively  large  amounts  of  a  carrier-free  radionuclide  of  Interest. 

The  assay  would  consist  only  of  measuring  a  well-defined  gamma 
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photopeak .  However,  it  would  still  be  necessary  to  show  for  each 
radionuclide  that  the  same  results  were  obtained  by  both  methods,  end 
this  would  again  necessitate  solving  the  above  problem  of  assaying  a 
small  amount  of  a  radionuclide  from  a  relatively  large  amount  of 
radioactivity.  In  addition,  it  would  be  desirable  to  assay  all  the 
radionuclides  of  interest  at  one  time  in  order  to  cut  down  the  total 
analysis  time  and  to  be  assured  that  all  the  radionuclides  underwent 
the  same  experimental  conditions  and  could  thus  be  directly  intercom- 
pared. 

A  direct  solution  of  this  problem  was  to  develop  a  computer  pro¬ 
gram  that  would  resolve  a  weakly  active  sub-mixture  of  gamma-emitting 
radionuclides  into  its  components.  There  are  several  requirements  and 
conditions  to  be  considered  for  any  computer  analysis  of  the  gamma 
spectra  of  a  sub-mixture  of  radionuclides:  (l)  All  the  major  radio¬ 
nuclide  components  present  must  be  identified,  (2)  the  total  amount  of 
activity  to  be  resolved  may  be  very  small,  (3)  the  inert  matrix  contain¬ 
ing  the  radionuclides  will  be  variable. 

There  are  various  computer  methods  presently  available  for  resolv¬ 
ing  gamma  pulse-height  distributions.  They  include  spectrum  stripping, 
solution  of  simultaneous  equations,  least  squares  fitting,  Iterative 
method?,  synthetic  methods,  and  linear  programming  (some  are  iterative). 
A  brief  outline  of  each  method  follows : 


(a)  Spectrum  stripping  means  Identifying  the  radionuclide  with  the 
highest  energy  peak  and  then  subtracting  its  normalized  pure  spectrum 
from  the  spectrum  of  the  mixture  and  repeating  this  procedure  until 
all  components  have  been  identified  and  quantitated. 

(b)  Simultaneous  equations  may  be  used  to  solve  a  spectrum  of 
identified  components  by  dividing  the  spectrum  into  regions  covering 
the  photopeaks  only.  Then  a  set  of  simultaneous  equations  is  developed 
relating  the  counts  observed  in  the  channels  in  the  photopeak  regions 
to  those  measured  in  a  set  of  standard  spectra  of  the  radionuclides. 

Within  the  various  channels  which  constitute  the  peak  the  differ¬ 
ent  radionuclides  at  "unit"  composition  are  known  to  contribute  ST^ 
counts  (i  Indexes  the  standard  radionuclide,  J  the  channel).  In  any 
one  channel  the  counts  in  the  mixture  spectrum  UHK(j)  are  given  by 

WK(J)  -  J  ST^ 


or  UKK(l)  -  a1  Sl^  +  a2  ST^  +  a3  ST^  + . »B  ST^ 

UWK(2)  -  «!  ST12  ♦  a2  STgg  + . •„  8 

a 

UlOC(m)  *  a,  ST.  +  . a„  ST 

x  ini  n  nm 


where  the  a1  are  the  unknown  composition  coefficients  of  the  various 
radionuclides  in  the  mixture.  These  sets  of  simultaneous  equations 
may  be  solved  by  straightforward,  if  tedious,  algebra. 
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(c)  Least  azures  fitting  extends  the  regions  considered  In  b 
above  Into  every  point  In  the  spectrum.  Thia  gives  as  many  equations 
as  channels;  each  equation  Involves  an  error  term  vhich  Is 

In  the  normal  vay.  That  la, 

m  n  2 

l  (»*<«-£•,«„  ) 

>1  1-1 

Is  minimized. 

(d)  Synthetic  methods  utilize  standard  spectra  to  give  Informa¬ 
tion  regarding  such  effects  as  backs  catter,  X-ray  production,  $  absorp¬ 
tion,  self -absorption,  summation  effects  and  detector  performance. 

for  a  given  energy  the  method  predicts  a  pulse-height  spectrum,  tm- 
approach  can  be  extended  to  analysis  of  a  composite  spectrum  by 
synthesizing  a  tried  composite  spectrum  comparing  It  with  the  experi¬ 
mental  spectrum,  and  then  modifying  the  synthetic  model  by  iterative 
procedures  until  a  match  Is  obtained. 

(e)  Linear  programming  utilizes  the  data  by  requiring  a  solution 
of  a  series  of  Inequalities  rather  than  ordinary  simultaneous  equa¬ 
tions.  One  method  of  solution  of  a  set  of  inequalities  (restrictive 
equations  containing  slack  variables)  is  called  the  simplex  method. 

For  example  instead  of 

•l  OT11  +  “2  ST2l  ‘  UNK(1) 
ve  have  ax  ST^  +  ag  ST^  s  UHX(l) 
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and  a  alack  variable)  r^,  la  Introduced  so  that 
a1  STU  +  a2  ST^  +  rx  -  UHK(l) 

These  canputer  methods  attempt  to  give  quantitative  guana  spec¬ 
trum  analysis.  All  the  techniques  are  sound  mathematically,  but  the 
Inherent  statistical  fluctuations  in  counting  lead  to  difficulties  in 
using  methods  a,  b  and  c  above.  Method  e  is  used  more  vhere  the  energy 
lines  of  gamma  emission  are  important  and  the  radionuclides  are  unknown. 

The  method  of  computer  analysis  developed  for  this  program  com¬ 
bines  some  of  the  features  of  a  and  e  in  order  to  meet  the  requirements 
mentioned  above  for  the  type  of  identification  needed  in  this  work. 


COMPUTER  PROGRAM  DESCRIPTION 

Briefly,  the  program  works  as  follows:  It  identifies  the  major 
peak  in  the  gamma  pulse-height  distribution,  compares  it  with  each  of 
the  major  peaks  in  a  library  of  pure  gamma -emitting  radionuclides 
under  the  same  geometry  conditions,  and  subtracts  a  small  fraction  of 
the  identified  pulse-height  distribution  from  the  unknown  mixture.  The 
major  peak  of  the  resulting  pulse -height  distribution  is  identified  and 
the  process  repeated  until  no  peaks  can  be  Identified  or  unaccounted- 
for  peaks  are  identified  which  are  not  in  the  library.  If  necessary 
these  unaccounted-for  pulse-height  distributions  are  put  into  the 


5 


library  ui  used  to  subtract  these  components  from  the  Mixture*  The 
residual  counts  in  the  pulse  height -distribution  are  analysed  by  a 
linear  programming  technique  (not  yet  used  in  present  experiment).  This 
method  is  particularly  suited  for  pulse-height  distributions  that  have 
feu  total  counts  (less  than  a  thousand)  or  have  poor  statistics  result¬ 
ing  from  the  stripping  program. 

Figure  1  shows  a  flow  diagram  of  the  computer  program.  Appendix  I 
®iT*a  the  program  in  detail.  This  program  represents  a  very  simple 
straightforward  analysis.  Many  refinements  can  be  added  to  it  as  re¬ 
quired  for  practical  applications. 

BCPEMKHTTAL 


Several  synthetic  mixtures  were  prepared  to  test  the  computer  pro¬ 
gram.  The  first  mixture  tried  was  composed  of  48.4  f  Ti^,  22.1  %  Cs13\ 
88 

•od  29*5  ^  Y  •  Each  ot  these  radionuclides  has  several  gamma  photo¬ 
peaks  over  a  2  Mev  range.  The  gamna  pulse-height  distributions  of  the 
separate  radionuclides  and  of  the  mixture  were  detected  with  a  3  in.  x 
3  in.  Hal(n)  crystal  and  analyzed  with  a  1024  channel  TMC  pulse-height 
analyzer.  The  computer  program  was  modified  to  plot  the  resulting 
pulse -height  distribution  after  50  1*  of  one  of  the  games  photopeaks 
was  subtracted.  Presently  the  major  peak  in  the  mixture  is  identified 
and  50  £  of  this  peak  and  a  proportional  amount  of  the  rest  of  the 

•These  are  larger  subtractions  than  would  normally  be  unde.  Normally 
only  a  few  percent  of  the  photopeak  would  be  subtracted  at  each  pass. 

The  plots  were  made  after  each  subtraction  in  order  tint  trouble  that 
developed  could  be  spotted  easily. 
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BEAD  BOW  BINARY 
DATA  CARDS 


Pig,  1  Flow  Diagram  for  Computer  Program  Used  to  Identify  Components 
in  a  Mixture  of  Gamma-Emitting  Radionuclides, 


pulse-height  distribution  of  ths  radionuclide  In  the  llbraxy  Is  subtrac¬ 
ted  without  reading  out.  mils  process  is  then  continued  until  no  more 
peaks  can  be  Identified  or  until  a  residual  pulse-height  distribution 
la  obtained  that  does  not  belong  to  the  library. 

RESULTS  AND  DISCUSSION 


The  results  of  a  test  of  resolving  a  mixture  of  the  three  radio¬ 
nuclides  Is  shown  In  Table  I  and  Figures  2-12.  Pulse-height  distribu¬ 
tions  of  the  pure  radionuclides  are  shown  In  Figures  13-15.  Each  pass 
represents  Identification  of  the  largest  peak  in  the  mixture  and  removal 
of  50  f  of  the  activity  in  that  peak  by  subtraction  of  an  appropriate 
amount  of  the  total  pulse-height  distribution  of  the  pure  Identified 
radionuclide  In  the  library.  Eleven  passes  were  made,  resulting  In  the 
removal  of  94  $  of  the  total  activity.  Machine  time  was  15  minutes. 

At  each  stage  the  remaining  peaks  could  easily  be  Identified  visu¬ 
ally  and  they  corresponded  with  peaks  in  the  library.  A  test  to  remove 
as  much  activity  as  possible  has  not  been  made  yet. 

The  small  errors  (0.5  -  2.0  percent)  in  the  analysis  are  considered 
very  encouraging  since  many  possible  refinements  to  the  computer  program 
have  not  been  made.  Two  of  the  refinements  are  concerned  with  peak 
identification  where  two  or  more  peaks  are  close  In  energy  and  with 
using  the  simplex  method  to  further  resolve  the  presently  remaining  re- 
eldue  of  the  gamma  pulse -height  distribution. 
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Table  1 


Stepwise  Subtractions  of  Radionuclides  Prom  a  Composite 

Mixture 


Pass 

Activity  Subtracted  at  Each  Pass  (cm) 

45 

Ti 

-73 

1 

66854 

2 

32927 

3 

35019 

4 

41703 

5 

20652 

6 

12974 

7 

9385 

8 

9799 

9 

4017 

10 

5489 

11 

1692 

Total 

counta 

subtracted  11487 5 

53482 

72354 

%  found 

47.7 

22.2 

30.1 

i  added 

48.4 

22.1 

29.5 

f  error 

-  1.4 

+  0.5 

+  2.0 
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It  appears  that  this  technique  offers  several  advantages  over  other 
compiter  methods  that  have  been  reported : 

(1)  The  high  and  low  energy  portions  of  the  pulse-height  distribu¬ 
tion  do  not  have  to  be  included  if  the  major  peaks  are  not  in  those 
areas.  The  low  energy  portion  is  very  sensitive  to  shifts  in  energy 
and  is  the  most  complicated  part  of  the  spectrum,  while  the  high  energy 
portion  often  contains  only  minor  fractions  of  the  gemma  emitters  (that 
must  be  included  in  the  usual  stripping  techniques).  Therefore,  the 
range  of  energies  can  be  quite  small,  if  desired,  and  the  detector  gain 
can  be  set  so  that  the  highest  energy  peak  will  determine  the  upper 
gain  setting  and  thus  maximize  the  resolution  of  the  detector  for  the 
gamma  photopeaks. 

(2)  This  computer  program  for  gamma-ray  analysis  is  very  sensitive 
to  pulse-height  distributions  with  very  small  numbers  of  counts.  The 
only  limiting  factor  in  the  program  is  whether  a  peak  can  be  identified. 
Peak  identifioation  is  determined  by  whether  the  total  counts  in  an 
energy  span  (corresponding  to  a  peak)  are  greater  than  the  total  counts 
in  any  other  energy  span  (that  might  be  due  to  perhaps  one  or  more  high 
statistical  fluctuations).  Even  large  statistical  fluctuations  can  be 
controlled  to  a  great  extent  by  arbitrarily  limiting  the  count  ratio  of 
neighboring  channels  to  less  than  a  certain  maximum  factor  such  as  would 
occur  at  the  inflection  point  on  the  high  energy  side  of  a  gamma  photo¬ 
peak.  Tests  are  being  made  for  mixtures  containing  less  than  103  total 
counts. 
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(3)  Although  It  is  required  to  have  in  the  library  all  the  radio¬ 
nuclides  in  the  mixture  to  get  maximum  aenaitivlty,  a  decided  advantage 
doea  exiat  in  being  able  to  identify  the  realdue  pulae-helght  distri- 
butlon  which  could  be  due  to  one  or  more  unexpected  radionudidea. 

With  accurate  knowledge  froet  the  literature  of  one  or  more  of  their 
gamma  photopeaks,  identification  ahould  be  eaay.  If  aueh  a  radionuclide 
la  easily  produced,  ita  apectrum  may  be  added  to  the  library,  and  if 

it  is  not  easily  produced  ita  pulse-height  distribution  can  be  accura¬ 
tely  synthesized  from  its  known  decay  scheme. 

(4)  Since  only  a  portion  of  the  pulse-height  distribution  of  the 
mixture  is  removed  at  each  pass,  the  true  percentage  of  each  radio¬ 
nuclide  present  ia  approached  asymptotically.  Thia  means  that  if 
components  are  not  present,  good  results  for  all  components  will  be 
obtained  when  a  large  fraction  of  the  residue  pulse-height  distribution 
still  remains.  In  the  example  given  above,  for  Instance,  when  only 
two  passes  have  been  made  for  each  radionuclide,  or  77  percent  of  the 

total  activity  removed,  the  error  for  the  Ti^  -  4.5  %,  Cs13^  -  15  $ 

88 

an“  ^  ■  7*5  When  minor  components  are  present ,  similar  results 

can  be  expected  when  two  passes  have  been  made  for  each  radionuclide. 

(5)  This  computer  program  should  do  well  for  radionuclides  having 
gamma  photopeaks  very  close  to  one  another  as,  for  example,  Zr^3  and 

95  ill!  lliji  1QQ  *1 

Nb  ,  Ce  and  Ce  ,  and  Ru  and  Ru  .  Advantage  can  be  taken  of 
the  fact  that  operations  are  being  performed  on  the  pulse -height 
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distribution  that  are  constantly  changing  its  shape.  For  example,  in 
an  eqjiilibrium  mixture  of  Zr95  and  Nb95,  the  peak  at  O.7U  Mev  might  be 
Identified  as  either  all  Zr95  (0.7a  Mev),  or  Nb95  (0.76  Mev).  Subtrac¬ 
tion  of  10  percent  of  the  peak  with  the  Zr95  staaiard  would  produce  a 
larger  shift  in  the  energy  of  the  resulting  peak  than  a  10  percent  sub¬ 
traction  of  Nb9*.  This  is  because  there  is  more  of  the  Nb95  photon  in 
the  equilibrium  mixture.  Therefore,  the  program,  after  trying  both 

subtractions,  would  recognize  this  smaller  shift  in  energy  and  the 

95 

Nb  would  be  subtracted.  The  remainder  of  the  peak  would  be  reduced 
in  a  similar  manner.  (The  50  percent  peak  subtractions  for  each  pass 
in  the  example  given  above  are  this  large  so  that  the  number  of  plots 
needed  in  the  development  of  the  program  would  not  be  too  large.  A 
much  smaller  percentage  subtraction  will  be  used  in  the  final  computer 
program.) 

(6)  Another  advantage  of  this  computer  program  is  that  very  com¬ 
plex  mixtures  of  radionuclides  do  not  make  the  analysis  much  more 
difficult.  As  long  as  one  photopeak  can  be  identified ,  the  other  ill- 
defined  portions  of  the  pulse-height  distribution  composed  of  other 
photopeaks  of  the  same  radionuclide,  ccmpton  edges,  etc.  will  become 
better  defined  since  the  components  containing  these  masking  details 
are  partially  removed.  For  example,  large  amounts  of  Te132-I132  in  a 
complex  mixture  could  completely  mask  a  small  amount  of  Cs137.  As  the 
Te13S-I132  was  removed,  the  Cs137  would  become  better  defined.  By  the 
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T*132  Ph0*0K"  “*  tta.  tot  of  the  C.W  lnter. 

"r1^  «“«*>»  I132  «»U  be  lee.  th»,  10  *rc.nt  of  the  c.»7 
photopeak . 

(7)  1M,  method  1.  not  very  ..nettle.  to  ,.i„  end  beeline  ehlfte 
of  the  detector.  The  only  llmt..,  fe.tor  u  tMt  ^  ;<  ^ 

larger  than  the  range  of  channela  f+  on 

20  in  fche  present  development) 

chosen  for  peak  identification. 

Many  refinements  may  be  made  in  +v,e  „ 

y  ce  made  in  this  computer  program  for  analya- 

in,  .more,  of  n*lonnelld.. .  ttl.  method  >m  ^ 

».lbl.  the  determination  of  „ry  l.r,,  of  l0.  „tivlty 

—dee  tbet  ere  .  nee..,..,  Mrt  of  the  e.ndy  of  t„ 

of  the  ch»le.l  „d  Phyelcal  .peed..  0f  «.r.et.r  „acl„  erpm.^ 
debris. 
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APPHTOIX  I 

Computer  Program  for  Identifying  Components  in  a  Mixture  of 
Oamaa-Qnltting  Radionuclides. 
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•  XEQ 

•  SYMTBL 

•  FORTRAN 

C  THIS  IS  A  STRIPPING  PROGRAM  FOR  GAMMA  SPECTRA 

DIMENSION  ST  f  9. 1032  )  .UNKI1032  )  ,  BUFFER! 1000) »  WTUNK  I  1032  I 
DIMENSION  STDI2000I 
NTHTRY  ■  0 
READ  100.NUMBST 
100  FORMAT  (111 

DO  3  I  ■  1.  NUMB ST 
CALL  REBIN  1 1024 .STD) 

DO  A  J  ■  1.1024 
4  STI I .J)  *  STD  I J) 

3  CONTINUE 

CALL  REBIN  I  1024. UNK) 

CALL  OREXEL  IUNX. NTHTRY) 

DO  9  LEE  ■  10.  1032 
XLEE  *  LEE 
FUDGE  «  2.0 

9  WTUNK ( LEE  I  «  SORTFIXLEE  I/FUDGE 
SIGUNK-0.0 

DO  10  LOVE  -  10.1032 

10  SIGUNK-S1GUNK  +  UNK  I  LOVE  I 

11  DO  12  JIM  -10.1032 
ARNIE  •  0.0 

12  ARNIE  -ARNIE  +  UNKIJ1M) 

IFIS1GUNK/ARNIE  -  25.1  13,13.14 

13  CALL  PEAKlt  JUNK.  UNK.  WTUNK I 

CALL  PEAKST  (JUNK.  LST .  IP,  ST,  UNK.  NUMBST.  NTHTRY) 

CALL  BRISS  (LST.  JUNK.  IP,  ST,  UNK) 

CALL  OREXEL  (  UNK, NTHTRY) 

GO  TO  11 

14  CALL  STOPLT 
CALL  EXIT 

END  10,1,0.0.01 


SUBROUTINE  REBIN  (NCHANL.  WD I 
DIMENSION  W0I1I.  ST0AI12),  ST0BI12I 
NCARD  -  (NCHANL  ♦  31)  /  12 
DO  50  IA  •  1,  NCARD 
IB  *  12  »  IA  ♦  1 
S  RCD 

DO  10  IC  *  1,  12 
S  CPY  STOAIICI 

S  CPY  STOBI IC) 

10  CONTINUE 
KOUNT  «  0 
DO  30  IC  «  1,  9 
K  »  10  -  IC 

S  CAL  STOBI ICI 

S  LRS  32 

S  T2E»30 

S  LRS  3 

S  TiE»21 

S  STO  STOR 

KOUNT  ■  K  »  1000 
S  CLM 

S  LDO  STOR 

S  21  LLS  1 

S  T2E*22 

S  STO  STOR 

KOUNT  «  KOUNT  ♦  K  «  100 
S  CLM 
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S  LOO  stor 
S  22  LIS  1 
S  TZ£*23 
S  STO  STOR 

KOUNT  «  KOUNT  *  K  •  10 
S  CLM 

S  LDQ  STOR 
S  23  LLS  1 
S  TZE*30 

FOUNT  =  KOUNT  +  K 

30  CONTINUE 

IF  (KOUNT)  60.  60.  31 

31  IF  (KOUNT  -  (A)  AO.  All  AO 

AO  PRINT  900.  (A.  KOUNT 

At  DO  SO  IC  *  1.  12 
10  «  II)  -  IC 
S  CLA  STOAI IC) 

S  FAD  STOA(IC) 

S  FSB  STOAI I C ) 

S  STO  WD( ID) 

50  CONTINUE 

60  IF  ( IA  -  NCARO )  61.  70.  70 

61  IB  »  NCARO  -  1 

DO  65  IC  •  IA.  IB 

DO  65  ID  =  1.  12 

IF  •  12  •  IC  +  ID  -  12 

65  WDIIE)  •  0.0 

70  RETURN 

900  FORMAT  ( 5H0CARD  15.  33H  IS  OUT  OF  ORDER.  CURRENT  CARD  IS  15.  19H 
1-  DATA  IS  ACCEPTED  ///I 
END  (O.U.O.O.O) 


•  FORTRAN 

SUBROUTINE  DREXEL  (UNK.NTHTRY) 

DIMENSION  SPI  1032I,  TITLE* 1*1 .  BUFFER  I  1000 ) .UNK ( 1032  ) 

IF  (NTHTRYI  111.111.112 

111  CALL  PLOTSiauFFERI 10001.1000) 

C  LOCATE  THE  PEN  AT  THE  ORIGIN  WITH  COORDINATES  (0.0. 0.0) 

112  CALL  PLOT!  5..  11..  -3  1 

CALL  PLOT!  0.0.  -10.375.  -3  I 

NOCHAN  =  102A 

NODCDS'A 

XSL-10.2A 

DC0L-2.A 

ORMEV  =0.0 

DXMEV-.l 

KEY  •  0 

INI  T  I A  -10 

NOSP  -1 

LAST  *102A 

3  CALL  PLOT  (-2.AA.  2. A.  3) 

CALL  PLOT  1-2. AA.  7.25.  2) 

C  THE  FOLLOWING  DO  LOOP  IS  TO  PRFVFNT  UNK  GETTING  CLOBBERED 
DO  2  1=1 . 102A 
2  SP  I  I  I  =  UNK  I  |  I 

C  DEFINITIONS  OF  ARGUMENTS  FOR  SUBROUTINES 
10  DELTA  =  XSL/  FLOATF (NOCHAN  ) 

XST AR T  =  FLOATF I  INITIA  -  1  )•  DELTA 
YSL  =  OCDL  *  FLOATF!  NODCDS  I 
YLEG  =  YSL  -  .9 
I F INAL  •  LAST  -  1NITI*  +  1 
00  50  N-l.NOSP 
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CALL  INPUT  (KEY. NOCHAN. SP.TITLEI 
IF  *  N  -  1  I  20.  20.  AO 
C  FIND  THE  MAXIMUM  VALUE  TO  BE  PLOTTED 
20  CALL  LARGCI  VMAX.  INITIA.  LAST.  SP  I 
C  SNOTN  SPLITS  A  NUMBER  INTO  SCIENTIFIC  NOTATION 
CALL  SNOTNI  VMAX .  S.  E  I 

C  1N0C0  IS  THE  EXPONENT  OF  THE  FIRST  DECADE 
INDCD  «  XINTF1EI  ♦  1  -  NODCDS 

C  GR I D  DRAWS  THE  SEMI-LOG  GRID  FOR  PLOTTING  THE  SPECTRA 
CALL  CRIDI  XSL. DCOL. NOCHAN. NOOCOS.  I NDCD.ORMEV.OXMEV  I 
C  SEMLGS  SCALES  THE  SPECTRA  TO  MATCH  THE  GRID 

*0  CALL  SEMLGSIO..DCi)L.lNOCD. INITIA, LAST.  SP  I 
C  HISTGM  PLOTS  THE  SPECTRA 

CALL  H1STGMI  XSTAR  T «  0.0,  DELTA.  1,  IF  INAL,  SP  I 
C  LEGEND  PRINTS  THE  IDENTIFICATION  OF  THE  SPECTRA  ON  THE  PLOT 
CALL  LEGEND! XSL+. 5.  YLEG,  TITLE  I 
50  YLEG  -YLEG  -1.25 

C  SET  UP  THE  ORIGIN  FOR  THE  NEXT  PLOT 

CALL  PLOT!  XSL  +  6,75.  0.0.  -5  I 
RETURN 

100  FORMAT!  '.15.  AF10.0  1 
ENDIO. 1.0, 0,01 


SUBROUTINE  PEAX1 (  JUNK ,UNK .WTUNK ) 
DIMENSION  UNKI1032) .WTUNKI1032I 
BIGUNK  •  UNKI 150  I 
DO  113  l«  151,1001 
IF  <1-1000)111, ni.m 

111  IFIBIGUNK-UNKI1I 1112.113.113 

112  BIGUNK  •  UNK! 1 1 
JUNK  ■  I 

113  CONTINUE 

114  RETURN 

END  10.1.0.0.01 


FORTRAN 

SUBROUTINE  PEAKST! JUNK.LST.I  .ST .UNK.NUMBST .  NTHTRYI 
DIMENSION  STI9. 10321 
NTHTRY  •  NTHTRY  +  1 
J  «  NUMBST  ♦  1 
IF  (  NTHTRY  -  60  I  5,80.80 
5  DO  60  I  *  1  ,  J 
IPASS  ■  0 

IF  I  I  -  NUMBST  I  10,10,80 
10  BIGST  •  ST  I  1,1501 
DO  15  K  ,  '51,  1001 

IF  (  BIGST  ST  I  I.  K  )  I  12.  15.  15 
12  BIGST  -  ST  I  I.KI 
LST  ■  K 

15  CONTINUE 

L2  ■  LST  ♦  50 
BIGST  -ST  «  I  .  L21 
DO  20  K2  ,  L2.1001 

IF  (  BIGST  -  ST  II.  K211  16,  20.  20 

16  BIGST  «  ST  (  I.  K2I 
MST  ■  K2 

20  CONTINUE 

CALL  OPCON  (2,6.51 

PRINT  82  .  ST!  I >2  I , LST. MST. JUNK 

82  FORMAT (9X.8HSTANOARD,F9»0,7HHI GH  PK.I9.6HLOW  PK.I9.9X.AHJUNK.  191 
24  DO  30  NP  ■  1,40 

IF  IJUNK-LST-20  ♦  NP)  30.40.30 


31 


1 


JO  CONTINUE 

IPA5S  »  IPASS  ♦ 

33  DO  33  JP  •  1 #40 

IF  ( JUNK-MST-20  ♦  JP  I  33.40*33 
33  CONTINUE 
60  CONTINUE 

AO  PRINT  41  .  ST( I .2  I .JUNK.  LST  *  MST 

41  FORMAT! 9X.2JHPEAKC0MPARI SON  FNO.ST* .F5.0. SHJUNK*. 1 3.4HLST* . 1 3. 1 5.  I 
IF  UPASS  1  70.70.69 

69  LST  •  MST 

70  RETURN 

60  PRINT  61.  NTHTRY 

61  FORMAT  I  9X  •  33HNO  PEAK  COMPARISON  FOUND  NTHTRY  «.1S) 

CALL  NOPCON 

CALL  STOPLT 

CALL  EXIT 

END  (  0.1.0.0.01 

SUBROUTINE  UKISS  (LST. JUNK. I  .ST.UNK1 
DIMENSION  UNK  (  10321.  S T ( 9 . 1032  I . TAKOPH ( 1032 ) 

19  »L$T-S 
JB  *LST+5 
SUMUNK*0.0 
SUMST  *0.0 

DO  143  KB  «  IB. JB 
SUMUNK  .SUMUNK  ♦  UNK I  KB  I 

145  SUMST  .  SUMST  +  ST ( 1  .KBI 
FACTOR  •  SUMUNK/ I  2 . *SUMST I 
DO  146  LB*  1 0 . 1032 

TAKOPH  (LB)  »FACTOR»ST(I  .LB) 

146  UNK I  LB  I  «  UNKILBI  -  TAKOPH! LB ) 

CONTINUE 

RETURN 

END  (0.1. O.O.O) 


SUBROUTINE  INPUTIKEYi  NOCHANfSP.  TITLE) 
DIMENSION  SPI  1032  |i  T|TL£<  16  I 
N  ■  NOCHAN  ♦  6 
40  TITLE C 1 )  •  SP ( 2 ) 

TITLE (2)  .  SP I  3 1 

TITLEI3I*  I960.  ♦  MODE  I  SPU)*.00l  «10.  ) 
TlTlEltla  .001*SP(S>  ♦  MOOFI  SP(4)  I  1000.  ) 
TItLEIS )»SP( l I  •  ,1 
IF (  key  I  SOI  601  70 
SO  TITLE 16)  ■  1000,*SP(6>  ♦  SPI?) 

DO  55  I  .  11236 
S3  SP(!>  ■  SPf  l  +  8  ) 

GO  TO  80 

60  T1TLEI6I  •  0.0 
00  63  I  ■  9IN0CHAN 
63  TITLE16)  •  T I T L E < 6 )  ♦  SPI  1 1 
GO  TO  80 

70  TITLE  <6 )  »  SPIKEY  ) 

60  CALL  SNOTNI  TITLEI6) •  T1TLEI7I.  TITLE  I  a )  I 
RETURN 

100  FORMAT (  14FS.0.  2X  ) 

ENolOil lOiOlO ) 
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SUBROUTINE  LARGE  component  finds  The  largest  component  of  A 


c  vector  bftween  the  i  and  l  components 
subroutine  largci  c*  i.  l.  v  i 
c  c  is  the  largest  component 

C  I  IS  THt  FIRST  COMPONENT  CHECKED 

c  l  is  the  last  component  checked 

c  v  is  the  name  of  the  vector 

dimension  v<ii 
c  «  vin 

DO  20  K  «  I.L 
IFI  V(K|  -  C  I  201  20»  10 
10  C  .  VIK) 

20  continue 

RETURN 

ENd(0»0»0»0*0  I 


SUBROUTINE  GRID  (XSL,  oCOLi  NOCHAN,  NODCDS.  1NDCD  ,ORMEV,  DXMEV  I 


c 


c 

c 

c 

c 

c 

c 


XSL  IS  THE  LENGTH  OF  THE  X  AXIS 

dcdl  is  the  length  of  a  decaoe 
nochan  is  the  number  of  channels 
nodcds  is  the  number  of  decades 
inocd  IS  the  exponent  of  the  value  of 


t  inches  i 

I  INCHE5  I 


The  v  ORIGIN 


draw  tear  line 

CALL  PLOT  (-2,**,  2,4,  31 
CALL  PLOT  (-2.4A,  7.25,  21 


DELTA  •  10. •  XSL  /  FLOATFI  NOCHAN  I 
YSL  *  DCDL  •  FLOATFI  NOOCDS  I 
IF(  1NDC0  +  NODCDS  I  5*5, 10 

5  5iLLSVMBUt",7J,<k‘3,,21,5H*  <11,  90, *5  ) 

GO  TO  J  5 

10  CALL  SYMBL*(-.7S«4.3,.21«5HA  III,  90., 5  | 


15 


LAX  IS  CRAWS 
CALL  LAXISI  0 


VERTICAL  LOG  SCALES  TO  FORM  SEMI-LOG  GRIDS 
. ,0, , XSL, DCDL, ,OSiNODCDS» INDCD  I 


PHAXIS  DRAWS  AND 
CALL  PHAXISI  0.0 


LABELS  THE  SCALE  ON  THE  X  AXIS 
,0,,  XSL  .DELTA,  .05  ,  ,0001  ,  JO,, 


CALL  SYMBLA  I  4«1«-.J,,21,10HCHANNEL»  I  ,0,,  10  ) 


‘1  ) 


-.05,  ORMEV,  DXMeVi 


call  PHAXISI  0.0,  YSL.  XSL. DELTA, 

IFI  DXMEV  I  20.  25,  20 

5»!'h  YSL+,!  ’  *21  *  12HENERGY  (MEVI  .  0.0. 

”  RETURN  PLAL+,S’  9,J5’*21,1jhU,s,n"r,0.L..  0.,l2) 
ENdIO, 0,0, 0,01 


2  I 
12  I 


33 


nrtrvnnonn 


SUBROUTINE  HISTGMI  X!  Y*  DELTA*  initia*  LAST*  Yh  ) 

x  is  the  x  coordinate  of  the  first  point  t 

Y  IS  THE  LOCATION  OF  THE  Y  COORDINATE  OF  THE  ORIGIN 

delta  is  the  distance  between  points  i 

INiTIA  IS  The  first  component  plotted 
last  is  the  last  component  plotted 
yh  is  the  vector 

dimension  YH(l) 

CALL  PLOT!  X*  Y*  3  I 
XHi  .  X 
C 

C  IF  DELTA  is  very  small  the  histogram  becomes  a  pt. 

1 F (  DELTA  -  .02  )  15*15*5 

5  DO  10  I  •  INITIA  *  LAST 
CALL  PLOT!  XHl.YH(I)  *  2  I 
C 

C  INCREMENT  THE  x  VALUE 

XHl  •  XHI  +  DELTA 
10  CALL  PLOT!  XH1*YH(I)  *  2  ) 

GO  TO  25 

15  DO  20  I  ■  INlTIA*LAST 

call  PLOT!  XH1*YhI||  •  2  I 

c 

c  increment  the  x  value 
20  XHi  ■  XhI  +  delta 
25  CALL  PLOT!  XHl*  Y.  2  I 
RETURN 

ENd  (0*0*0.0*01 


LAXIS  DRAWS  A  VERTICAL  LOG  AXIS 


C  X  AND  Y  ARE  THE  COORDINATES  OF  THE  OF  THE  GRID  ( 

c  xsl  is  the  length  of  the  x  axis  ( 

c  dcol  is  the  lingTh  of  a  decade  I 

c  t jci_  is  the  length  of  the  Tic  marxs  ( 

c  nodcds  is  The  number  of  decades  desired 
c  indcd  is  the  exponent  of  the  initial  decade 

subroutine  laxisix*y.xsl«ocdl*ticl*nodcds*indcd  i 

dimension  dyiioi 

XE  »  X  -  *63 

YE  .  Y  -  .11 

XEE  «  X-  *21 

YEE  ■  Y  ♦  .06 

INDEX  ■  INDCD 

XX  .  X-  TICL 

XSL  ■  X+XSL 

XXX  •  XSL+  TICL 

L  .  NODCDS  +  1 

c 

c  compute  vertical  increments  within  decades 

DO  5  J*  2*10 
3  ovtj)  .  DCDL  «  LOGlOFf  FLOaTF(J)  1 
CALL  PL0TIX.Y.3I 

C  draw  left  axis  AND  DECADE  LINES 

Yd  «  Y 

DO  15  J  *  1 .NODCDS 
DO  10  K«  2*  10 
YY  ■  YD  *  DYI  X) 

CALL  PLOT  I  X.  YY*2> 


INCHES  ) 

I  INCHES  ) 
INCHES  ) 


TO  PT.  PLOT 


INCHES  ) 
INCHES  I 
INCHES  I 
INCHES  ) 


c 


c 


CALL  PLOT  (  XX.YY.2I 
10  CALL  PLOT  (X.  YY*2) 

CALL  PLOT  ( XSL iYY* 2 ) 

CALL  PLOT  (X.  YY.2I 
13  YO  •  YD  +  OCOL 

DRAW  VERTICAL  labels 
DO  20  J  »  It  L 

call  NUMBER  (XE  »YE  •• 21*10. *0.*  -1) 
rx  •  FLOATF!  INDEX  I 
CALL  NUMBER  (XFE*YEE*.07«FX  .0,1  -1) 
INDEX  »  INOEX  ♦  1 
YE  s  YF  ♦  DCDL 
20  Y-E*  YEE*  DCDL 

draw  right  axis 

CALL  PLOT  (XSL.  Y  *  3  I 
YD  «  Y 

DO  30  J  »  ltNODCDS 
DO  25  K«  2*10 
YY  *  YD  +  0Y(  X) 

CALL  PLOT  (XSL.  YY l  21 
CALL  PLOT  (XXX.  YY i  2) 

23  CALL  PLOT  (XSL.  YY.  2) 

30  YD  •  YD  ♦  DCDL 

CALL  PLOT  (  X.Y.3I 
return 

end  (0.0.0.0.01 


SUBROUTINE  LEGEND!  X.  YLEG*  TITLE  ) 
DIMENSION  TITLEI16I 
Y  .  YLEG 
XP  »  X  ♦  1. 

XT  •  X  ♦  2.75 

GC  •  TITLE (6l 


CALL 

SYMBLA 

I  X.Y. 

•  14* 

1 0HSAMPLE 

CODE. 

0.. 

18  ) 

CALL 

number 

(XP.Y. 

•  u* 

T I TLE I 1 1 .  0..  -1 

1 

call 

number 

(XT.Y. 

•  14* 

T I TLE (2  I  *  0,.  -i 

1 

Y.  I 

.25 

CALL 

symbl* 

(X  .Y. 

•  14* 

17HYEAR 

DAY. 

0.. 

17  ) 

CALL 

NUMBER 

(XP.Y. 

•  14* 

T I TLE ! J I •  0.*  -1 

1 

CALL 

number 

(XT.Y. 

•  14* 

TITLE  1 A  1  *  0,.  3 

1 

Y. 

v-  .23 

CALL 

SYMBLA 

(  X.Y. 

•  14* 

17HC0UNT ING  INTERVAL. 

0,. 

IT  1 

CALL 

number 

(XT.Y. 

•  14* 

T jTLEf J)»  0,.  1 

1 

Y  . 

Y  -.25 

CALL 

SYMBLA 

(  X.Y. 

•  14* 

26HGR0SS  COUNT 

X10 

•0. *26) 

CALL 

NUMBER 

IX4*1.73#Y», 

1A.T[TLE<T)»  o,.  6 

CALL 

NUMBER! 

XU, 73* 

,07.  .07»T|TLE($I> 

o«*  -1 

) 

return 

end  (0.0. 0.0. 01 


pretty  horizontal  axis 


C  X  and  Y  ARE  THE  COORDINATES  OF  THE  ORIGIN  I  INCHES  ) 

c  xsl  is  the  length  of  the  axis  <  inches  i 
c  delta  is  the  distance  between  the  tic  marks  (  inches  i 

c  tic  is  the  lfnc.th  of  the  short  tic  marks  (  inches  i 
c  orchan  is  the  base  value  for  computing  the  labels 

c  ox  is  added  to  orchan  at  every  long  tic  and  the  sum  drawn 

C  N  is  THE  N  USED  IN  SUBROUTINE  NUMBER  OE  THE  CAL  COMP  SERIES 

C 
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SUBROUTINE  PhAXISI  X.  Y.  XSL.  delta*  Ttc*  ORCMAN.  OX*  N  ) 

XA  •  X 

chan  ■  orchan 

YT  •  V  -  TIC 

DT  •  TIC  ♦  TIC 

aTic  »  ABSF  (  OT  ) 

nolab  »  XINTfI  ,2*XSL/OELTa  ♦  .001  ) 

CL*b  •  1.3  *  ATic 
YDT  ■  Y  -  DT 
YLab  ■  TOT  -  DT 
c 

C  START  pen  at  ORIGIN 

call  plot  I  XA  i  Y  *  3  ) 

c 

DO  20  x  ■  1  •  NOLAB 

DO  lo  J  •  1  .  * 

XA  ;  XA  ♦  DELTA 

c 

c  make  FOUR  SHORT  TICS 

CALL  PLOT  I  XA  .  Y  •  2  ) 

CALL  PLOT  (  XA  *  YT  *  2  ) 

10  CALL  PLOT  (  XA  .  Y  *  ?  ) 

C 

c  hake  and  label  long  tics 

XA  .  XA  ♦  DELTA 

chan  ■  CHAN  +  ox 

XLaB  *  XA  -  CLAB 

CALL  PLOT  I  XA  .  Y  *  2  I 

CALL  PLOT  (  XA  *  YDT  *2  I 

IF (  DX  I  IS. 20. IS 

15  CALL  NUMBER  (  XLAB  .  YLAB  .ATIC.  CHAN  »  0.0  .  N  ) 

20  CALL  PLOT  I  XA  .  Y  •  J  ) 

CALL  PLOT  I  XSL  .  Y  l  21 
C 

RETURN 

end  (0. 0.0.0. 0) 

SUBROUTINE  SEMLGS  (  Y.  DCOL.  INDCD.  INV.  LASTV.V) 

C  SEMI-LOG  SCALE  SHIFTS  .  LOGS.  AND  SCALES  A  VECTOR  FOR  SEMI-LOG  PLO 
C  Y  IS  THE  LOCATION  OF  THE  ORIGIN  I  INCHES  ) 

C  DCOL  IS  the  LFNGTH  of  A  DECADE  I  INCHES  ) 

c  indco  is  the  exponent  of  the  value  of  the  origin 

C  INv  IS  THE  INITIAL  COMPONENT  SCALED 

c  lastv  is  the  last  component  scaled 

C  V  IS  THE  VECTOR  TO  BE  SCALED 

c 

dimension  vui 
c 

VMJN  ■  10.»*INDCD 

BASE  ■  DCOL  »  ELOATFl  INDCD  I  -  Y 

J  •  1 

SF  .  .A3*29*5»DC0L 

c 

c 

DO  20  K  «  INV.  LASTV 
C 

C  the  if  STATEMENT  AVOIDS  UNDEFINED  OR  UNPLOTTABLE  VALUES 
I F  (  V  I X  |  -  VMIN  I  10.  10. IS 
10  VUI  •  Y 
00  TO  20 
C 

15  VUI  •  SF*  LOGFI  V < K I  I  -  BASE 
C 

20  J  *  J+l 
return 

ENoIO. 0.0. 0.01 
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this  subroutine  will  find  the  magnitude  of  a  number 


GIVEN  X  THE  ROUTINE  FINDS  S  AND  E  SUCH  THAT 

subroutine  snotni  x*  s»  e  i 

J  g  0 

S  .  ABSFI  X  I 

IFI  S  I  6*  6*  1 

1  IFI  10,  -  S  |  2.  2*  3 
!  5  i  S*  il 

J  a  J  ♦  1 

GO  TO  1 

3  IFI  S  -  1,  |  4,  i,  5 

4  S  .  S  *  10. 

J  .  J  -  1 

GO  TO  3 

5  S  .  S  IGNf I  S.X  ) 

6  E  ■  FLOATFI  J  | 

RETURN 

end  IOiO.O.O.O) 


•  S»10.«*E 
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UNCLASSIFIED 
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